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Abstract 


The compatibility conditions for the strain measures of the theory 
of micromorphic elastic solids are presented. Specific forms are 
given for the linear theory and for the linear micropolar elasticity. 
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1. Introduction 


Er ingen and Suhubi [1, 2] gave a theory of microelasticity in which 
the local intrinsic motion of material elements is taken into account 
in the continuum behavior of elastic solids. This theory and its 
later extensions and specializations show that a large variety of 
physical phenomena, hitherto not encountered in the classical theories 
of elasticity and fluid dynamics, can be treated. For example, the 
effects of couple stress, distributed body couples, surface tensions, 
and distributed microcracks can be formulated mathematically and the 
underlying physical phenomena can be studied systematically. Other 
work related to this has appeared both in the United States and Europe. 
The aim of the present paper is not an assessment of existing literature, 
but to supply the compatibility conditions missing in these works. 

For a discussion of various theories, we refer the reader to [3] and 
[4]. 


2 . Strain Measures 

According to the theory of micromorphic materials^, the motion of a 

material volume element AV is described in terms of the "average 

motions" of large numbers of microvolume elements AV , (a = 1, 2,...) 

contained in AV. The position vector of a material point in AV^ 0 ^ 

( a ) __ (ct) 

is marked by X = X + = where X is the position vector of the 
~ ~ _ (a) ~ 

center of mass P of AV and = is the relative position vector of a 

(ct) ~ 

material point of AV from the center of mass P, Fig. 2.1. The 
motion of the assemblage contained in AV is postulated to be 

x (a) = x(X,t) + X K (X,t) £ ( “ ) (2.1) 
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This terminology was introduced by Eringen [5] in his later work. 
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where repeated indices indicate summation over the range (1, 2, 3). 

In (2.1) x(X,t) is the spatial position p of the center of mass P 

at time t and the remaining terms on the right of (2.1) represent 

(a) 

the intrinsic motion of the elements AV about the center of mass. 

It is assumed that both "macro" and "micro"-motions 

x = x(X,t) , = x K (X,t)H^“' ) (2.2) 

are continuous and possess continuous partial derivatives with 
respect to X^ (the rectangular coordinates of X) and time t as many 
times as we desire. It is also assumed that the inverse motions to 
(2.2) exist in some neighborhood of X and are given by 

X = X(x,t) , 5 (a) =^ k (x,t)C^ a) (2.3) 


except possibly at some singular surfaces, lines and points. 


This is secured by assuming that the jacobians not vanish 

3x 

J = det (~) + 0 


j = det ( XkK > + 0 


(2.4) 


in some neighborhood of X during the time interval under consideration. 

Here det means determinant and x k £ are the components of referred to the 
spatial rectangular coordinates x^. Similarly the components of^^ referred 
to the material rectangular coordinates X^ are denoted by i.e.. 


*K = X kK ( ^’ t) ^k 




(2.5) 
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where i^, (k=l, 2, 3) and (K=l, 2, 3) are the unit base vectors 
of the spatial and material corrdinates x^ and X^. respectively. 
Deformation gradients (x, , X, ) and the microdeformations (X.., 9C ^ ) 

play important roles in the theory. Among them we have the 
orthogonality relations 


’‘k.A.Jl 6 k£ ’ 


X K,k X k,L 6 KL 


( 2 . 6 ) 


X kK®k£ 6 k2. 


^Kk X kL 6 KL 


(2.7) 


of which (2.6) is an immediate consequence of the chain rule of 
differentiation and (2.7) follows from carrying one of (2.2)^ and 
(2. 3) 2 into the other. 

By considering one set of the deformation gradients and one set of 
the microdeformations unknown, we can solve for the other sets from 
the nine linear equations (2.6) and (2.7). Thus for example 


cofactor x. „ 

y = kjJC 

J 


_ 1 _ 

2J £ KLM e k£m X £,L X m,M 


cofactor x 


& 


kX 


Kk 


1 

2 j e KLM e k«.m X «.L X mM 


( 2 . 8 ) 


(2.9) 


where e^ ^ and e^ are the alternating symbols. For the jacobians 
we also have 


J 


1 

6 S KLM e k5>m X k,K X £,L X m,M 


( 2 . 10 ) 


. = 1 

3 6 e KLM S k£m X kK X £L X mM 


( 2 . 11 ) 
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By differentiation from (2.10) and (2.11) we also have the identities 


3 J 


3x 


k,K 


cofactor of x, „ = J X., . 

k,K TC.k 


( 2 . 12 ) 


= cofactor of (2.13) 

Of these, (2.12) is the well-known identity of Jacobi while (2.13) 
is appropriate to the present theory. 

The differential elements in the undeformed and the deformed bodies 
are calculated respectively by 

» 

dX(a) = dX + dH = (X, k +# £jk ^)dx k +3t k de k (2.14) 

dx(a) = dx + d£ = (x, K + x L + XjrdH^ (2.15) 


where we drop superscripts a from H and and write X(a) and 
(a) (a) K k 

x(a) for X and x for simplicity. In this form X(a) and x(a) 
may be understood to denote the position vector of a material point 
in AV^ 01 ^ and in Av^ a ^ , respectively, when the number of such microelements 
in AV and Av is dense. 

The square of arc length ds(a) in the deformed body can be calculated by 



-1 


ds 2 (a) (C^ + + rj^Fg^Cj.g — N ) 


+ 2 ^KL + ' i, SL I ’RMK C RS"M' >dX K d:: L 


+ ' 1, ML >i, NK C MN d “K d “L 


(2.17) 


where 


C KL ( ?’ t) = X k,K X k,L 

(2.18) 

= X k,K X kL 

(2.19) 

r KLM < '?’ t) = X k,K X kL,M 

(2.20) 


-1 

the deformation tensors of this theory and 
given by 

C_, T is the inverse of 

IxL 

-1 


C KL = ^.A.k 

(2.21) 


-1 

Here C and C are the Green and Piola deformation tensors of classical 
continuum mechanics and i and F are microdeformation tensors of the 
theory of micromorphic materials. In this theory, therefore, to calculate 
length and angle changes, we need two sets of new strain measures 
namely V and r in addition to C of classical continuum mechanics. 


Upon introducing the displacement vector U and the microdisplacement 

tensors 0 by 
KL 


X k,K (6 LK + U L,K )6 kL 


^kK 


= ( 6 . 


LK 


+ 


6 kL 


( 2 . 22 ) 
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The deformation tensors (2.18) to (2.20) can be written as 


C KL = 2E KL + 6 KL 6 KL + U K,L + U L,K + U M,R U M,L 


5 ^KL + 6 KL = 6 KL + + U L,K + U M,K $ ML 


(2.23) 

(2.24) 


r = $ + u $ 

KLM KL,M N,K NL,M 


(2.25) 


where E and C VT , so defined, along with T are the material strain 
KL K.L K.LM 

measures of this theory. 


The entire theory can be developed, equally well, in terms of the spatial 
deformations and strain measures (cf. [2]). 


We note that there exist 6 + 9 + 27 = 42 strain components which are 
expressed through (2.23) to (2.25) in terms of 3 + 9 = 12 displacement 
components and . Clearly then any arbitrary set of C, ? and T do 
not correspond to physically realizable displacement fields U and $ . 

More precisely we may pose the question: What are the conditions for 

which a set of C, f and T give rise to single-valued displacement fields? 
This question is answered in the following section. 
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Equations of Compatibility 


The compatibility conditions are the restrictions that must be placed 
upon C, ¥ and T so that single-valued continuous displacement fields 
U and $ exist satisfying (2.23) to (2.25). In other words they are the 
integrability conditions for the set of partial differential equations 
(2.23) to (2.25). In principle, these conditions can be obtained by 
elimination of the gradients of U and $ from these equations by successive 
differentiation and elimination. This method is, however, hopelessly 
awkward. As an alternative method we make use a theorem of Riemann 
namely: For & symmetric tensor a^ _to be_ a_ metric tensor for 

a_ e uclidean space , it is necessary and sufficient that a^. be^ a_ non- 
singular , positive definite tensor and that the R iemann-Christof f el 

(a) 

tensor R, „ formed from it vanish identically. We have 
k£mn z - 


(a) 1 

R k£mn ” 2 ^ a kn,£m + a £m,kn a km,£n a £n,km^ 


-1 

+ a rs ([£m,s] [kn,r ] 

-1 

where a^ is the inverse matrix to a^ 
of the second kind defined by 


• [£n,s] [km,r]) (3.1) 

and [k£,m] is the Christoff el symbol 


[k£ ,m] 


2 (a km,£ + a £m,k & k£,m ) 

-1 -1 

a a„=a a„ = 6 
ns s l ns si ni 


(3.2) 


To apply this theorem we need to recognize the metric tensor in 

the deformed and undeformed bodies. Clearly Y and T cannot be 

KL KLM 

used for this purpose for in general V is not a symmetric tensor and 
of course T is not a second order tensor. We observe, however, that 
(2.17) may be expressed in another form namely 

ds (a) = ^ B dX A dX B , (A,B =1, 2,..., 6) (3.3) 
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r 


where we set dH. = dX. , dH„ = dX n , d5„ = dX- and 
1 4 2 5 3 6 


-1 


^*KL ^LK C KL + 2r (KL) < ' a) + r MK (a) r Nl/°^ °MN 
CcZ> _ Ct? _ 

^KK' ~ ^K'K ~ \ K'-3 + K’-3 F RK C RS 


r. . =r 


E T ¥ 

'v* ‘ot ' 1 T 


K'L' L'K' MK'-3 N L’-3 MN 


r KL (a) F KMl“m, 


(3.4) 


(K,L,M,N = 1, 2, 3), 


(K ' , L', M\ N' = 4, 5, 6) 


where a parenthesis enclosing indices indicates the symmetric part 
of the quantity involved, e.g., 


r (KL) 2 ^KL + r LK ) 


l - 3 > i! 


In (3.3 ) c? is a nonsingular symmetric positive definite second order 
AB 

tensor of the euclidean six-dimensional space. Therefore, for 

AB 

to be a metric tensor of this space, according to a theorem of Riemann, 
we must have 


R.ABCD = ° (A ’ c » D » = !» 6 > 

Substituting in place of a in (3.1) we get 

1 <€r + . if ) 

2 v AD,BC BC ,AD AC ,BD BD,AC ; 


(3.5) 


(3.6) 


-1 


( [BC,F] [AD,E]-[BD,F] [AC,E]) = 0 
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where 


[BC.F] 5 | (^ Br>c + 


& - 'S’ ) 

W CF,B BC ,F' 



(A,B,C, . . . ,F = 1,2,.. .,6) 


(3.7) 


These results are valid for curvilinear coordinates. 

The tensor can have a total of 21 independent components. In six 

dimensions the total number of nonvanishing components of R 

AJdCD 

is 105. Thus (3.6) gives a total of 105 compatibility equations for 
the deformation tensor These equations are second order partial 

differential equations. By use of the definitions one can express 
these equations in terms of C, ¥ and T. The resulting expressions are 
too complicated to reproduce here. Instead we give the explicit 
forms of the compatibility conditions for the linear theory. In this 
case all the nonlinear terms in (3.6) are dropped reducing it to 


^AD,BC 


+ £ 


BC ,AD 


^AC,BD 


BD , AC 


= 0 


(3.8) 


For the linear theory from (3.4) we have 



where 


(3.9) 


^KL ^LK “ 6 KL + 2E KL + Zr (KL) (a) 

^KK’ = ^K'K = 6 K K ' -3 + K K'-3 + r K'-3 K (a) 

^K'L' = ^L'K' = 6 K'-3 L ' -3 ~ 2E K'-3 L'-3 + 2 <? (K'-3 L'~3) 

where E, £ and T are the material strain measures of the linear 
theory defined by 
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(3.11) 


E KL S I <\,l. + °L,K> 


= $ + u 

KL KL L,K 


r = $ 

KLM KL,M, 


r KL r KML £ M 


The symbols 6 , _ and 8 , , are the usual Kronecker deltas which 

K. K. — j K — j L — j 

assume the value 1 when the indices take the same numerical values and 
zero otherwise, e.g., 


6 K K'~3 * 


1 when K'-3 = K 


0 otherwise 


We note the following identities for future use 

k KL,K' = r K K'-3 L + r L K'-3 K > SfKL.K'L’ = ° 
^KK', L' = P K'-3 L ' -3 K ^KK'.L'M' = ° 

^Tk'l' ,m’ = 0 


(3.12) 


Now decompose (3.8) into the following set of independent equations 


e* c** _ 

^KL,MN + 6 MN,KL “ ^KN ,ML “ ^ML.KN “ U 

^KL,MK' + ^K' ,KL ~ ^KK' ,LM ” ^LM,KK' = ° 

^KL.K'L' + ' L ' , KL “ ^KL'.LK' “ ^LK',KL' = ° 

^KK'jL'M' + ^L'M' ,KK' " ^KM'.K'L' “ ^K'L',KM’ “ ° 
^K'L'.M'N' + ^M'N'.K'L' " ^K'N',L'M' " ^L'M'.K'N' = ° 


(3.13) 
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Upon substituting from (3.10) and using the identities (3.12) we find 
that the last two equations of (3.13) are satisfied identically and the 
first three reduce to 


e kl,mn + e mn,kl e kn,lm e lm,kn 

4- i f fr +■ r 'i + (r +r 

2 U KSL 1 LSK' ’MN v MSN NSM' ’KL 

(r KSN + r NSK ):, ML " ^MSL + ^LSM^ ’KN^ “S 

£- - P + r + r 

°MN,KL C 'KN,LM KNL,M LNK,M 

” r LNM,K " r MNL,K + ^NSMjKL “ F NSK,LM^"s 
2 £(MN),KL ~ 2E MN,KL " r MNL,K ~ r NMK,L = ° 


(3.14) 


(3.15) 


(3.16) 


Hence equations (3.14) to (3.16) are the compatibility conditions of 
the linear theory of micromorphic elastic solids. When E is not 
considered to be arbitrary we have eighteen independent equations 
instead of 105. When Eg is considered arbitrary then the coefficients 
of Eg must be set equal to zero independently. This leads to 


e kl,mn + e mn,kl e kn,lm - e lm,kn ° 

(v + r 'i 4-Cr + r } 

v KSL LSK' ’MN v MSN NSM' ’KL 

_ (r + r 'i _ rr +r =o 

U KSN ISr’ML U MSL LSM' ’KN (3.17) 

^MN,KL _< ^KN,LM + F KNL ,M + F LNK,M “ F LNM,K “ F MNL,K = ° 

r NSM,KL ~ r NSK,LM = ° 

o > pp — r — r = n 

C (MN) ,KL MN ,KL MNL ,K NMK ,L 
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Upon using (3.17)^ we find that O.l?)^ is satisfied identically. 
The integration of (3.17)^ gives 


r NSM,K ~ P NSK,M ^NSKM 


(3.18) 


where ft Mt!MK . is independent of X and H and is subject to the 


NSMK 
restrictions 


fi NS (KM) ° 


(3.19) 


With the help of (3.18), equations (3.17)^ and (3.17),. can be 
written as 


( ^MN,K ^KN,M + r KNM " F MNK ) ’L + °KNML + ^MNLK + ^LNMK ° (3 - 20) 


2[< ^(MN),K " E MN,K ~ r (MN)K ], L °MNKL 


(3.21) 


These are integrated to give 

SmN,K _( ->KN,M + r KNM ~ r MNK = “ ( \nML + ^WlK + ^LNMK^^L + B KMN 


(3.22) 


' (MN) ,K E MN,K " r (MN)K 2 ^MNKL ^ + 


KMN 


(3.23) 


where B,„ _ T and D T „„ T are independent of X and - and are subject to 
KMN KMN r ~ J 


B (KM)N °* 


d kmn d knm 


(3.24) 


If the displacement field U,. possesses continuous second order partial 

derivatives and the microdisplacement field $ possesses continuous first 

KL 

order partial derivatives then 


^KLMN B KLM °KLM ° 


(3.25) 


be seen by merely substituting E , E and r from (3.11) 

K.L KL KLM 


as can 

on the left of (3.22) and (3.23). 


13 



Hence 


e kl,mn + e mn ,kl e kn,lm _ e lm,kn 

F -F + r - r =o 

C’MN.K °KN,M KNM MNK 
^ (MN) ,K " ” r (MN)K = ° 

r - r =o 

KLM , N KLN ,M 


(3.26) 

(3.27) 

(3.28) 

(3.29) 


and the proof of 

Theorem 1 : The necessary and sufficient conditions for the integrability 

of the system (3 . 11) for _a simply-connected domain is the satisfaction 
of the compatibility conditions (3 . 26) t<D (3 . 29) . 

When the compatibility conditions (3.26) to (3.29) are satisfied then 

single-valued solutions of the system (3.11) exist for U and $ and 

K. KL 

have the most general form 

U K = U K + Vl + \ 

$ KL = $ KL + ^KL 


(3.30) 

(3.31) 


where and B R are independent of and and 

fi (KL) = ° (3 ' 32) 

Physically + B^ represents a rigid body displacement field. 

The compatibility conditions corresponding to the eulerian counterpart 
( \£’ £ ki and Y k£m } ° f the Straln raeasures E KL ’ ^KL and r KLM haV6 


14 



identical forms to (3.26) and (3.29) where capital letters are replaced 
by small letters and r is replaced by -y. For the definitions of the 
eulerian strain measures and other related accounts the reader is referred 
to [1], 


4 . Compatib i lity Conditions for Micropolar Elasticity 

For micropolar elasticity we have,' by definition 

$ = — <£) 

KL LK 

This suggests introducing an axial vector $ such that 

K. 

$ = — -p d) = — p cf> 

K 2 KLM ML ’ KL KLM M 


(4.1) 


(4.2) 


The linear parts of the strain measures given by (2.23) to (2.25) now 
read 


E 

<5 

r 


KL 

KL 

KLM 


2 (U K,L + U L,K^ 


U - e $ 
L,K KLM M 


— e $ 

KLN N,M 


(4.3) 


From these equations it is clear that 



J KL 


r (KL)M 


(4.4) 


Using these we see that (3.28) is satisfied identically. Thus: 

Theorem 2 : The necessary and sufficient conditions for the integrability 

of the system (4.3) for a_ simply-connected domain is the satisfaction of 
the equations (3.26), (3.27) and (3.29). These are the compatibility 
conditions of the linear theory of micropolar elasticity, Eringen [6]. 
(3.26) and (3.27) were obtained by Sandru [7]. 
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